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Abstract 

We prove, using the coordinate Bethe ansatz, the exact solvability of a model of three 
particles whose point-like interactions are determined by the root system of 32 ■ The statistics 
of the wavefunction are left unspecified. Using the properties of the Weyl group, we are also 
able to find Bethe equations. It is notable that the method relies on a certain generalized 
version of the well-known Yang-Baxter equation. A particular class of non-trivial solutions 
to this equation emerges naturally. 
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Using the coordinate Bethe ansatz [T], C.N.Yang [2] solved a model of n particles, with un- 
specihed statistics, interacting via contact interactions. This procedure led him to discover the 
celebrated Yang-Baxter equation, found also, in the context of statistical physics, by R.J.Baxter 
[3]. The potential used in the approach of C.N.Yang is intimately linked to the simple root system 
of sin- The generalization to the other root systems has also been intensively studied (see for 
example 111 El) In the case of the so„ and sp„ root systems, a new type of equation, the so-called 
reflection equation mm, is obtained. It plays a fundamental role in the study of integrable system 
with boundaries. 

In this letter, we demonstrate, using the procedure of [2], the exact solvability of a model 
based on the 02 root system. This gives rise to a generalized version of the Yang-Baxter equation 
niEiini which is peculiar to 02 amongst the simple Lie algebras and is distinct from the usual 
Yang-Baxter and reflection equations. We then hnd the Bethe equations for the model, which 
involves some interesting subtleties. 

The Model 

Let us consider a system of three particles with positions Xi, X 2 and X 3 , whose interactions are 
specihed by the Hamiltonian 

^ ^ q2 3 ^ ^ ^ ^ 

^ = ~^'^ + ‘^9s^5{xi-Xj)+2gL - -{xi + xj)). ( 1 ) 

i=l * *,i=l i,j,k=l 

Yi i¥=j¥=k^i 

Here gs, gx are real parameters characterizing the strength of two types of interactions. Physically, 
the gs term is the usual contact term - specifying how particles interact when they collide - while 
the gi term may be thought of as describing a contact interaction between each particle and the 
centre of mass of the remaining pair. We wish to solve the spectral problem 


H(j){xi, X 2 , X3) = E(t){xi,X 2 , X3). (2) 

The motivation for the Hamiltonian (CD is that it is related to the root system of the exceptional 
Lie algebra 02 [5], just as the model of n particles with purely contact interactions is related to 
that of sin- the root system of 02 is 


A = {cj - Cj, Ci + Cj - 2efc|l ^ i ^ j ^ k ^ i ^3}, 


(3) 


where {cj} an orthonormal basis of In terms of these roots, the Hamiltonian reads as 


H 



+ ^9{a,a) S{{a,x)) 


( 4 ) 


where x = 92 = QSi 9& = 9l and (■, ■) is the usual scalar product. 
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Coordinate Bethe Ansatz 


To translate the problem into one which may be solved by the coordinate Bethe ansatz, it is nsefnl 
to consider the Weyl gronp associated to the Lie algebra 92 |1]. To each root a G A is associated 
a reflection in the hyperplane Tia perpendicnlar to a\ 

s„(/?)=/?-2 (5) 

(a, a) 

The set of these reflections, together with the identity /, form the Weyl gronp of 92 , which maps 
the set of roots to itself. The gronp is generated by the reflections in the two simple roots 

«! = ei - 62 and 0^2 = 62 + 63 - 2 ei, ( 6 ) 

which we denote by, respectively, T and R. The relations obeyed by these generators 

{Tf = I , {Rf = I and {TRf = / , (7) 

completely specify the gronp, which is isomorphic to the dihedral gronp Dq. 



Fignre 1: Root system of 92 , showing the labelling of Weyl chambers. 


The twelve hyperplanes T^q, define twelve domains, called Weyl chambers, in These cham¬ 
bers are characterized by 

: 0 < Xs^l - Xs^2 < - Xs^l , ( 8 ) 
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where Xsj, is the 
equivalent to 


component of the vector (sa) ^(^r). Let us remark that conditions 


: Xso < Xs^i and x, i < 


Xsc» + 


are 

(9) 


In each Weyl chamber, the three particles are in a given order, as are the centre of mass of the 
two extreme particles and the middle particle. The importance of the dehnition of the Weyl 
group in this context is two-fold |1]: the hyperplanes, Tia, are the domains of the conhguration 
space where the interactions take place; and the action successively of R and T allows us to 
describe, starting from one region Ws„, the eleven other regions. (See hgure QJ) Using these 
two properties, the spectral problem m may be written equivalently as a free Hamiltonian, for 
{Xi,X2,X3) ^ {Ha I a e A}, 


-E 


■^(j){Xi,X2,X3) = E(l){xi,X2,X3) 


( 10 ) 


and boundary conditions on hyperplanes Ha (« ^ A) 


4^\{a,x)=0+ 4^\{a,x)=0~ and {c^i'^)(p\{a,x)=0+ ((®) E‘2g(^a,a)') 4^\{a,x)=0~ (H) 

where V = equivalence between the (5-potentials in the Hamiltonian Q and these 

boundary conditions (dU is well-known (see for example 0). 

We now determine the eigenfunctions of the Hamiltonian by solving equation m with bound¬ 
ary conditions (dH). We make the following ansatz for 0: in the region Wq (with Q G Dq), the 
eigenfunction of the Hamiltonian is written as follows 


X2, X3) = ^ exp{i{kp,XQ))Ap{Q). 
PeDe 


( 12 ) 


This ansatz is similar to the one suggested by H.Bethe 0: the only difference is that here the 
sum is over a dihedral, rather than permutation, group. These eigenfunctions obviously satisfy 
relation dni) with 

(13) 

i=l 

We need to determine the parameters Ap{Q) present in the ansatz so that boundary conditions 
m are satished. Because of dehnition (jH)) of the region Wq {Q G Dq), the eigenfunction 0 q 
adjoins two boundaries: Xqi = and Xqi = xq 2 . The boundary conditions imply the 

following constraints between the vectors Ap (whose components are Ap{Q), P,Q E Dq): 

Apr = B{{kp,a 2 ))Ap and App = Y{{kp,ai))Ap (14) 


where 


B{k) 


kR + igp 
k - igi 


and Y{k) 


kf + igs 
k - igs 


(15) 
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and the operators R and T, which provide a realization of the gronp Dq, are dehned by 

RAp{Q) = Ap{QR) and fAp{Q) = Ap{QT) . (16) 

One can now nse relations m to calcnlate recnrsively all the Ap starting from, for example, 
Aj. However, dne to the relations m satished by the generators of Dq, the following consistency 
relations appear between the operators m 

B{k) = , Y{k) = {Y{-k))-^ (17) 

and 

y(fci - k2)B{2k^ -k2- k^)Y{k^ - k^)B{k^ + k2- 2ks)Y{k2 - k^)B{2k2 -k^- k^) 

= B{2k2 - fci - k^)Y{k2 - fc3)5(fci + k2- 2 fc 3 )H(fci - k^)B{2k^ - k2 - k^)Y{ki - ^ 2 ) (18) 

Using solely that R and T satisfy relations dlD, one can verify by direct compntation that these 
equations hold, hnishing our argument about the exact solvability of the 02 model. The relations 
m are the usual unitarity relations whereas (HHD is a generalization of the Yang-Baxter equation 
niHiini. We discuss the implications of its appearance briefly below. 

Bet he Equations 

In order to hnd the Bethe equations we suppose that the three particles live on a circle of hnite 
circumference 2L. One must specify carefully how the particles interact in this case. As on the 
infinite line, there are interactions due to direct collisions of particles. But, since on a circle there 
is no preferred notion of which particle lies on the left, on the right, or in the middle, it is now most 
natural to assume that each of the three particles interacts with the midpoint of the remaining 
pair. What is more, for each pair there are really two “midpoints”, and we assume that the third 
particle interacts with both. (The midpoint of two particles would otherwise jump discontinuously 
as they pass opposite points, which seems physically unappealing.) 

To make contact with the preceding section we need a prescription to identify the set of possible 
conhgurations of the particles on a circle with the group Dq. We make the following choice. We 
“cut” the circle at the point opposite the mid-point of the closest pair of particles (call them A 
and H - of course, the particles cannot be equally spaced for then each would lie on the mid-point 
of the other two). We then unwrap the circle to get an interval of length 2L, and the positions of 
the particles on this interval correspond to a unique element of Q G Dq, just as in the previous 
section. (Figure El illustrates an example in which Q = I.) 

By construction, neither A nor B can he at an endpoint of this interval (in fact, neither can 
be closer than 2L/3 to an endpoint) so we need only specify the boundary conditions for the 
remaining particle, C.^ On examining the definitions of the various regions, one sees that C can 

^Note also that the prescription ensures that - to take the example in figureEl- 1 reaches the “internal” midpoint 
(23) before 2 has a chance to reach the “external” midpoint (13). So the boundary condition on particle C (3, in 
this case) is genuinely the only new condition, not present in the model on the line. 
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(12) (21) (23) (12) 



Figure 2: A configuration in region W/. (The midpoint line (13) is not shown.) 


reach only the right endpoint for all conhgurations corresponding to 

Q G := {/, T, RTR, TRTR, RTRT, TRTRT} 


( 19 ) 


and that it can reach only the left endpoint for all conhgurations corresponding to 

Q e B+W = [RTRTR, {RTf, TR, R, TRT, RT} (20) 

where we dehne W = RTRTR. More specihcally, when one starts in a region Wq, with Q G 
and moves particle C through the right boundary, one reaches the region Wqw- As it crosses the 
boundary, C interacts with the “opposite” (on the circle) midpoint of A and B, so the conditions 
are continuity with a jump in the hrst derivative, just in m above. Explicitly, one hnds 


4^\{ao,XQ)=-2L- <P\{ao,XQ)=+2L- 


( 21 ) 


and 

{oiQyWQ')(j)\(^ao,XQ)=—2L~ ((® 0 ) R 4^\{ao,XQ)=+2L~ 1 ( 2 ^) 

where ao = ei + €2 — 2e3 is the lowest root. Note that W = RTRTR is the rehection in the plane 
orthogonal to ciq, so that 

{ao^xq) =-{ao.XQw) (23) 

and hence one can replace Q with QW in these boundary conditions without altering their content 
- as expected, since this is supposed to describe the boundary between Wq and Wqw- In fact 
the boundary conditions may be re-written most symmetrically as 


011 


ao,XQ)=-2L- 


I {oiO,XQy^r) = — 2L 


( 24 ) 
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and 


Q)(l^\{ao,XQ)=-2L - q) +‘^9l) (l)\{ao,XQw)=-2L-■ (25) 

Substituting the Bethe ansatz for 0 (in region Wq on the left and region Wqw on the right), 
one has, from the hrst boundary condition, 

PS-De PS-De 

Now the vectors n = :^(1,1,1), «i and are orthogonal, and hence 


{kp,XQ) = {kp,n){n,XQ) + ^{kp,ai){ai,XQ) + ^{kp,ao){ao,XQ); 
furthermore W hxes n and cki and inverts ckq. Thus 

^i{kp ,n){n,XQ)+^{kp ,ai){ai,XQ) 

P&B+ 


(27) 


(28) 


(29) 


= ^ gi(fcp,n)(n,XQ) + §(fcp,Qi)(ai,XQ> + Apw{QW)e^'~^^’"’°'°'‘ 

PeB+ 

which can hold for general x only if, for all P G , 

7lp(g)e-f + 24pvK(g)e+^<^^’“°> = 7lp(giV)e-^<^^’“°> + 7lpv^(giT)e+^<^^’“°\ 

Meanwhile, from the boundary condition on the hrst derivative of (j) one hnds, by very similar 
reasoning, that 

t{kp, ao)Ap{Q)e-f^>^^’^°^ - t{kp, ao)Apw{Q)e+f^^^’^^^ 

= {-t{kp, «o) + 2gp) Ap{QW)e-f^'^^’^°^ + {t{kp, «o) + ^gp) Apw{QW)e+f^’^^’-°'> (30) 

(here we used {kp,XQw) = {kpw,XQ) for the hrst term on the right). 

After eliminating Ap-\^[QW) in the second equation using the hrst, one arrives at 

( 31 ) 

{kp, ao) — igp 

where hhAp(g) = Ap(ghT). This is true for all P G but there is some redundancy, for if 
(pTT|l holds for P it also holds for PT : 

e^^kPTPo) = e^^kPTPo)j^p^^ 

= e~^^^’°‘°'^Y{{kpw,0'i))Apw 


\^(fcp,ao)bT + ic/p ^ 

= T((/cpiy,ai))—-^^-Ap 

\kpi olq) — igp 

{kp,aQ)W Pigp 


{kp, ckq) — igi 
{kpT,ao)W Pigp 
{kpT, cxo) — igL 


Y{{kpw, 0!i))Ap 
Apt 


(32) 


^In fact it is, as an immediate consequence, true for all P G Dq: on setting P = P'W one finds the same 
equation for P'. 
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(the essential point is that oq and ai are orthogonal, so hh = hxes ai, T = Sa^ hxes oq, and 
[hh, T] = 0). It therefore suffices to consider 

P e {/, RTR, TRTR}. (33) 

Using now the relations (d, which come from the “interior” boundary conditions, m yields 
three equations for Aj{Q). For example, in the case P = I, 

Aw = Artrtr = B{{kRTRT, o: 2 })y {{kRTR, ai))B{{kRT, (^ 2 ))Y {{ku, ai))B{{k, a 2 ))Ai. (34) 


To make the content of the resulting equations clearer, it is helpful to dehne some new operators. 
First, for every root a, let SaAp{Q) = Ap{Qsa) and dehne 


Za{k) 


{k) Cr) T i(j(a,a)^a 

(/c, Cr) '^9(a,a) 


(35) 


which has the property that Za{k) ^ = Za{—k). This notation is compact but rather opaque, so 
it is useful to dehne also 


5 ' i 2 — -^£ 1 - 62 ) ^23 — Ze2-e3, S^i — ( 36 ) 

(whose inverses we write as S 21 , S 32 , and S' 13 ) and similarly for the long roots 

- 2 £ 1 , K3I = ^£ 3 +£ 1 - 2 £ 2 ) -^12 = ^£ 1 +£ 2 - 2 £ 3 - ( 37 ) 

In terms of these operators we hud, for P = RTR, P = TRTR and P = I respectively, 

= ^lAi := S2iik)KUk)-^K^23ik)KUk)-^S3iik)K^,,ik)Ai, ( 38 ) 


^23 = 


g^(fc,£3+£i-2£2>^^ = ^271/ := Kl,{k)-^KUk)KUk)-^S32{k)KUk)Su{k)Ai, (39) 

^2^(fc,,,+,,_2£3>^^ = ^371/ := K^,,{k)-^S,3{k)KUk)S23{k)KUk)-^KUk)Ar. (40) 

These are the Bethe equations for the problem, and the task is to show that the commute. 
But, before this, it is important to observe that there is further redundancy. The exponentials on 
the left hand sides contain the three long roots of 02 , which are of course co-planar and indeed sum 
to zero. Thus there are really only two independent equations, and, multiplying the equations 
above together, we hnd that we must have 


'^l'^2'^3^/ — A-I- (41) 

To see that this is in fact true - which is a good consistency check - and to verify the commu¬ 
tation relations = 0, one use the following properties of S and K: 

= Kf^Su (42) 
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Kl, {Kl,)-' Kl, = Kl, (A'?3)-' A'fj 


( 43 ) 


SijA-'jSia (A-fj) ■' SjaA-fj = A'JjSaa (Af^)SiaA'jSia. (44) 

These may be verified directly. (See also P [Till ITT| .l 

The reason there are only two equations - even though there are three momenta ki - is that 
we chose to apply periodic boundary conditions in a fashion which made no reference to any 
fixed point on the circle. (It is more usual [121121 ESI to take, for example 0|a;i=o = 0 Ui= 2 L, but 
in the present case the nature of the interactions make this technically inconvenient.) Thus the 
symmetry of the problem under rotations was kept manifest throughout and the corresponding 
conserved quantity, the total (angular) momentum 

P = ki + k 2 + k^ = \/3(fc, n), (45) 

dropped out of the calculation. But of course, as in any quantum-mechanical system of particles 
on a circle, P is the generator of rigid rotations, and the invariance of the problem under a rotation 
through one complete turn produces the quantization condition 

Y = = g2iL(A:,ei-|-e2+e3} (46) 

This, together with any two of ()d8ll4()j) . gives the complete set of quantization conditions on the 
momenta ki. 

It is interesting to note that, at least in the centre of momentum frame ki + k 2 + k^ = 0, the 
equations (IdHHdOj) have the intuitive interpretation one would expect: for example, (Ei becomes 

^2^Lk,A^ = ir2yA;)-iir2^(A:)i^32(A:)-'^32(A;)i^?3(^)5i2(fc)kl/ (47) 

and describes the process of moving particle 2 clockwise through one complete revolution while 
the other particles remain fixed. Thus, the first event is the scattering of 1 and 2 (hence S' 12 ) 
followed by particle 2 interacting with the midpoint of 1 and 3 (giving scattering with 3 

( 532 ), and so on. 

Conclusion 

To conclude, let us comment briefly on the generalized Yang-Baxter equation (I44j) we obtained. 
Like the Yang-Baxter and reflection equations, (m may be represented diagrammatically. This 
is shown in Figure 121 where, to simplify the picture, we restrict ourselves to the case where 
Xi+X 2 +X‘i = 0. The three arrows represent the “trajectories” of the three particles and the double¬ 
line that of their centre of mass. The intersection between the two arrows a and b corresponds to 
the scattering of Sab of particles a and b. The intersection between the arrow a with the double- 
line corresponds to the scattering ( 6 , c ^ a) between a and the centre of mass. Obviously, we 
recover the usual representation of the Yang-Baxter equation by removing the double-line. This 
occurs in the limit Ql = 0 , for then = 1 . 
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Figure 3: Pictorial representation of the generalized Yang-Baxter equation iH) 

Since the Yang-Baxter and reflection equations play a fundamental role in the development of 
integrable models and quantum groups, it is natural to speculate that the generalized Yang-Baxter 
equation dH might also have interesting applications. In particular, we hope that they will allow 
one to study integrable models where the interactions between three particles are not factorisable. 
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